Tidal tails around globular clusters. 
Are they a good tracer of cluster orbits? 
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ABSTRACT 

We present the results of detailed N-body simulations of clusters moving in 
a realistic Milky Way (MW) potential. The strong interaction with the bulge 
and the disk of the Galaxy leads to the formation of tidal tails, emanating from 
opposite sides of the cluster. Some characteristic features in the morphology 
and orientation of these streams are recognized and intepreted. The tails have 
a complex morphology, in particular when the cluster approaches its apogalacti- 
con, showing multiple "arms" in remarkable similarity to the structures observed 
around NGC 288 and Willman 1. 
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Actually, the tails are generally good tracers of the cluster path quite far from 
the cluster center (> 7-8 tidal radii), while on the smaller scale they are mainly 
pointing in the direction of the Galaxy center. In particular, the orientation of 
the inner part of the tails is highly correlated to the cluster orbital phase and 
to the local orbital angular acceleration. This implies that, in general, the or- 
bital path cannot be estimated directly from the orientation of the tails, unless 
a sufficient large field around the cluster is available. 

Subject headings: globular clusters: general, galaxies: kinematics and dynamics, 
methods: n-body simulations 



Introduction 



It is commonly accepted that the present globular cluster (GC) population in the Milky 
Way represents th e survivor of an initially m or e numerous one, dep auperated by many dis- 
ruptive processes jMurali & Weinberg! [l997al jbl: iFall fc Zhangll200lh . 

The observational results obtained in the last decade have clearly confirmed the role played 
by the Milky Way environment on the GC evolution. Up to now, about 30 galactic GCs 
show evidences of star depletion due to the tides caused by the field of the Galaxy. 



The fi rst evidence of the existence of tails surrounding GCs were achieved by lGrillmair et al. 



( 119951 ). Using colour- magnitude selected star counts in a dozen of galactic GCs, these authors 
showed that in the outer parts the stellar surface density profiles exceeded significantly the 
prediction of King models, extending also outside the tidal radius. They identified this sur- 
rounding material as made up of stars in the act of being tidally stripped from the cluster field 
and pointed out the importance of defining the large-scale distribution of extra-tidal stars on 
the sky to obtain constraints and traces for the GC orbits. More recently , other works con- 
firmed and enlarged Grillmair et al.'s findings (jLehmann fc Scholzl 119971 ; iTesta et al.l |2000| ; 
Leon et al.ll2000l : ISiegel et al.ll200ll : iLee et al.ll2003r) givi n g evid ence of the existence of many 
GCs surrounded by haloes or tails. For M92, ILee et ah] (120031 ) found, also, that the extrati- 
dal material is not randomly distributed, being the density profile in the tails shallower for 
bright stars than for fainter ones. 

This was the state of the art until the spectacular findings of two tidal tails emanating 
from the outer part of the Palomar 5 globular cluster and covering an arc of 10 degrees 
on the sky, corresponding to a (proj ected) length of 4 kpc at the distance of the cluster 
( lOdenkirchen et al.ll200ll . |2002| . 120031 ) . obtained in the framework of the Sloan Digital Sky 
Survey (see also http://www.sdss.org). The stellar mass in the tails of this sparse, low-mass, 
halo cluster (with an estimated concentration parameter c = 0.7) adds up to 1.2 times the 
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mass of stars i n the cluster, estimated in t he range between 4.5 x 1O 3 M and 6 x 1O 3 M . 
More recently, iGrillmair &: Dionatod (120061 ). still using SDSS data, have detected a contin- 
uation of Pal 5's trailing tidal stream out to almost 19 degrees from the cluster. Combining 
this with the already known southern tail of Pal 5 yields a stream some 9 kpc long on the 
sky. 

Substantial tidal streams have recently been found associated with another low-mass and 
low-concentration GC in the SDSS area: NGC 5466. For t his cluster iBelokurov et al.l (20061) 
reported giant tails extended for about 1 kpc in length and Grillmair fc Johnson (12006 ) . still 
using SDSS data, suggested that a « 13 kpc tidal stream of stars, extending from Bootes to 
Ursa Major, could also be associated with thi s system. 

Toge t her with the so-called Sagittarius stream (IMateo et al.lll998t lYanny et al. hood IMartmez-Delgado et 



200ll : iMajewski et al.. 1 1200.4 IMartmez-Delgado et alihooi T" which emerges from a dwarf 



galaxy that is currently being accreted by the Milky Way, Palomar 5 and NGC 5466 rep- 
resent outstanding examples of ongoing tidal erosion of stellar systems in the Milky Way, 
being, up to now, also the only two globulars known for which such extended stream-like 
structures have been detected in the Galactic halo. 



One of the first numerical invest igations of the role played by a galactic tidal field on 
spherical stellar systems was that of iKeenan &: Innanenl (119751 ). who studied the effect of 
realistic, time-varying tidal fields on the stellar orbits in a star cluster. They numerically 
integrated the equations of motion of three-bodies in models of spherically symmetric clus- 
ters which, in turn, move in eccentric orbits in the field of a model galaxy. One of the main 
conclusions of this work, which extended previous investigations made by iKingl (119621 ). was 
that star clusters rotating in a retrograde sense are more stable in a tidal field than clusters 
with either direct rotation or no rotation due to the contribution of the Coriolis acceleration, 

acting in the same direction as the gravitational attraction for retrograde motion. 

More recent works on weak t idal encounters based on a Fokker-Planck approach (lOh &: Lin 
1992 : Lee &: Goodman 1995 ) and self-consistent N-body techniques ( Grillmair 19981 ) con- 
firmed that the interaction with external tidal field, combined with two-body relaxation in 
the core of the cluster and following replenishment of stars near the tidal radius, causes a 
flow of stars away from the cluster. The stripped stars remain in the vicinity of the cluster 
for several orbital periods, migrating either ahead of the cluster or falling behind, giving rise 
to a slow growth of the tidal tails. 

A semi-analytic study of the development of tidal streams in galactic satellites was done by 



Johnston! ((1998), who gave estimates for the rate of growth of tidal tails. The effects of a 



realistic galactic t idal field (including b oth bulge, halo and disk) on GCs were investigated 
few years later by I Combes et al.l (119991 ). The main findings of the work were the following: 
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stars escaped from the system go to populate two giant tidal tails along the cluster orbit; 
these tails present substructures, or clumps, attributed to strong shocks suffered by the clus- 
ter, and are prefe rentially formed by low mass stars. 

Yim fc Led (120021 ) performed N-body simulations of GCs orbiting in a two-componen t galaxy 
mode l (with bulge and halo and no disk), using the direct-summation NBODY6 code (lAarseth 
19991 ) and focusing their attention, in particular, to the correlation between tidal tail elonga- 
tion (described by mean of a 'position angle', defined as the angle between the direction of 
the tail and the galactic center direction) and the cluster orbit. They found that, on circular 
orbits, tidal tails mantain an almost constant position angle (~ 60°), while GCs on non 
circular orbits show a variation of the position angle, according to orbital path and phase. 
The position angle increases when the cluster heads for perigalacticon. On the other hand, 
it tends to decreases when the cluster heads for apogalacticon. 

Finally, some authors investigated also the dynamical evolution of some globular clusters 
in the tidal field of the Galaxy. In this context iDehnen et al.l (120041 ) modelled the disruption 
of the globular cluster Pal 5 by galactic tides. Pal 5 is remarkable not only for its extended 
and massive tidal tails, but also for its very low mass and velocity dispersion. In order to 
understand these extreme properties, they performed many simulations aiming at reproduc- 
ing the Pal 5 evolution along its orbit across the Milky Way. They explained the very large 
size of Pal 5 as the result of an expansion following the heating induced by the last strong 
disk shock about 150 Myr ago. The clumpy substructures detected in the tidal tails of Pal 
5 are not reproduced in their simulations, so that they argued that these overdensities have 
been probably caused by interaction with Galactic substructures, such as giant molecular 
clouds, spiral arms, and dark matter clumps, which were not condsidered in their modeling. 
These simulations also predict the destruction of Pal 5 at its next disk crossing in about 110 
Myr, suggesting that many more similar systems have once populated the inner parts of the 
Milky Way but have been transformed into debris streams by the Galactic tidal field. 
In this context, it may be interesting mentioning the recent numeric al work devoted to the 
study of smaller size systems (open clusters) in the MW tidal field ( IChumakfe Rastorguev 
20061 ). which confirms already known results in the case of the external part of GC tidal tails) 
abou t the alignement of stars of the tidal stream around a common o rbit in the external 
field jGrillmairill998l : ICombes et alill999l : ICapuzzo Dolcetta et al.ll2005h . 



In the above sketched theoretical and observational background, this work -which 
i s inserted in a wider study on the dynamics of globular clusters i n exte r nal tidal fields 



(ICapuzzo Dolcetta et all l2005t iDi Matteo et all l2005t Miocchi et al.l l2006al : Miocchi et al. 
2006bh - is devoted to clarifying the connection among tidal tails and cluster orbit. We will 



describe and discuss the mechanisms that determines the tails morphology and how they 
depend on cluster trajectory and orbital phase. 
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For this purpose, we performed detailed iV— body simulations (with N = 10 5 ) of GCs moving 
in a realistic three-components (bulge, disk and halo) Milky Way potential. 
The paper is organized as follows: in SectfS] the Galaxy (Sect JUT] ) and cluster (Sect 12.21) 
models adopted are presented, so as the numerical approach used; in SectE] we deal with 
the main results of our work, showing the formation and development of tidal tails around 
the cluster (Sect. 13. ip . giving a qualitative approach for describing the tail morphology 
(Sect 13. 2p . presenting the numerical procedure adopted to fit tails direction (Sect J3.3P and, 
finally, discussing the tail-orbit alignement and its dependence on the orbital phase. In 
SectJU all the results are summarized and discussed. 



Models and methods 



2.1. Galaxy model 



The model adopted for the Galactic mass distribution is that of lAUen Sz Santillanl (119911 ). 
It consists of a three comp onent system: a spherical central bulge and a flattened disk, both of 
Miyamoto fc Nagail (119751 ) form, plus a massive spherical halo. The gravitational potential is 
time indipendent, axysimmetric and given in an analytical form which is continuous together 
with its spatial derivatives. 

Choosing a reference frame where the (x,y) plane coincides with the MW equatorial plane, 
the three components of the potential have, in cylindrical coordinates, the form: 
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where square brackets indicate the difference of the function evaluated at the 100 kpc and 
the generic r = \/ R 2 + z 2 extremes. The parameters in the formulas above are listed in 
Tabled 
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2.2. Cluster model and numerical method 



The total initial mass of the cluster was chosen Mqc — 4.7 x 10 5 M , i.e. a value com- 
pati ble with mass es of galactic globular c luster s lying inside 4.5 kpc from the MW center 



Harris ( 1996 ) and Pryor Sz Meylan ( 19931 )). The stellar mass spectrum of our clus- 



(see 

ter initial model was chosen as a Kroupa IMF sampled in the range O.1M — 2OM and 
" evolv ed" (in the sens e of ac c ounting for mass loss on the base of stellar evolution, accord- 
ing to IStraniero et al.l (119971 ); iDominguez et al.l (119991 )) up to 3 Gyr (which corresponds to 
our assumed cluster age). In this i nterval of time all masses grea ter than 1.2 M go into 
the O.5M — 1.2M range. As in ICapuzzo Dolcetta et al.l (120051 ). we sampled this mass 
function into 12 mass classes, equally spaced in a linear scale, whose space a nd velocity 



distribution were evaluated by the adoption of a multimass King distribution (jKinglll966 
Da Costa &: Freeman! 1 19 761 ). 

Obviously, the choice of N = 10 5 as total number of cluster stars, together with the given 
value of Mqc and of the above described mass function implied a rescaling of the star masses. 
Finally, to investigate the role played by the degree of cluster concentration, we considered 
clusters with two different values for the King concentration parameter c (listed in Table [2]). 



The clusters move on the y — z coordinate plane (the (x, y) plane corresponds to the 
Galactic disk), along orbits of different eccentricity, defined as 



r a + r p 



(5) 



being r p and r a , respectively, the GC pericentric and apocentric orbital distances. See Table 
[3] for GC orbital parameters and FigfTj for a plot of the different simulated orbits. 

All the simulations were performe d by means of the TreeATD co de, developed by two o f 
us ( IMiocchi &: Capuzzo-Dolcettall2002l ) which, as the original code by lBarnes fc Hutl (119861 ). 
relies upon a tree-algorithm for the gravitational force evaluation on the large scale and 
on a direct summation on the small scale. The code was parallelized to run on high per- 
forma nce computers (via MPI routines) e mploying an original parallelization approach (see 
again (IMiocchi fc Capuzzo-Dolcettall2002l . see)). The time-integration of the 'particles' tra- 
jectories is performed by a leap-fro g algorithm that uses individual and va riable time-steps ac 



cording to the block-time scheme (lAarsethlll985l ; lHernquist fc Katzlll9 89). See IMiocchi et al. 



(j2006al ) for further details. The only constant of motion that can be significantly used to 
check the quality of the orbital time-integration is the total energy of the cluster, i.e. in- 
cluding the contribution of the external field to the GC potential energy. We saw that the 
upper bound of the relative error in the energy conservation, AE/E, is 10~ 4 over the whole 
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simulation time. This is no more than one order of magnitude worse than the error we got 
in a comparison simulation for the same GC in abscence of external field. 



3. Results 

3.1. Formation and evolution of tidal tails 

In Figs. I2"f51 and HI the formation and subsequent evolution of tidal tails around the GC 
is shown, for all the simulations performed in the case of the least concentrated of the two 
GC models (the one with c = 0.81 in Table [2]). In the case of the more concentrated initial 
model cluster, the development of the tidal tails shows the same time dependence being just 
slightly less populated. 

After about 5 x 10 7 yr, tidal tails have clearly formed. They continuously accrete by the stars 
leaving the cluster, so that after 1.6 x 10 8 yr, in the case of the quasi-circular orbit (orbit 
III), they are elongated for more than 1.5 kpc each. In general, the extention of the two 
cluster tails depends on the velocity of the cluster and its variation along the orbit. For the 
quasi-circular orbit (orbit III) this velocity is almost constant and the two tails have nearly 
the same length. For more eccentric orbits (orbit I and II), the tail extending between the 
GC and the pericenter is more elongated than the opposite tail, because the stars in the 
latter have smaller velocities than those belonging to the former one. In any case, the tail 
that precedes the cluster always extends slightly below the orbit, whereas the trailing one 
lies slightly above the orbit in agreement with what was observed for Pal 5. This feature can 
be explained considering the role played by the Coriolis acceleration, as we will see later. 

As regards the tails orientation, these figures (Figs [2|3|4|) clearly show the alignment of 
the outer part of the tails (at distances from the GC center > 7-8 tidal radii) with the orbit. 
On the contrary, the alignment of the inner part is strongly correlated to both the orbit 
eccentricity and the GC location. Notice, indeed, that the inner tails are roughly aligned 
with the orbital path only when the cluster is near the perigalacticon of the more eccentric 
orbits (I and II, respect i vely s hown in Fig f5] and Figj3]). This confirms previous results in 



Capuzzo Dolcetta et al.l (120051 ). where analogous features were found for GCs movings in a 



triaxial potential. 

It is also worth noting the peculiar morphology in the streams when the cluster is 
approaching the apocenter in the eccentric orbits (orbits I and II): tidal tails divide into 
"multiple arms" , two nearly elongated along the GC path, while the other two point to the 
galactic center (an d anticenter) direction. At this regard, the tails observed around NGC288 



( ILeon et al.ll2000l ) have a similar complex morphology, showing three different "arms". In- 
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terestingly, this cluster has been estimated to be near its apogalacticon (IDinescu et al.lll999l ) 
in agreemen t with our results. A lso the galactic satellite Willman 1 shows multi-directional 
stellar tails (jWillman et al.l 120061 ). which according to our results, could represent evidence 
that this system is approaching the apocenter on an eccentric orbit. 



3.2. Tidal tails elongation: a qualitative description 



As anticipated in |Pi Matteo et al.l (120051 ). the simplified scheme of the formation and 
shape of tidal tails around GCs can be understood by the motion of a star escaping from 
a cluster moving in a spherical potential. Let us consider a rotating reference frame with 
the origin in the cluster center-of-mass, the (x',y r ) plane coinciding with the orbital plane 
and x' pointing to the galactic center. The equation of motion of the iy, star belonging to 
the cluster that moves around the galaxy center with variable angular velocity u> is (see also 
Appendix IA1 and FigjSJ) : 



r- = fj - r GC -2wxr[-o;x(wxr|)-wx r-, 



(6) 



where r' { is the position vector in the rotating frame, while r^ and the position vector of 
the GC center-of-mass, tgc-, refer to the inertial frame. For a star escaping through one of 
the unstable Lagrangian points along the x' axis, the first, second and fourth terms on the 
right-hand side of Eq. [H]are directed along the x' direction, while the third term (the Coriolis 
acceleration) and the fifth term are along y'. The latter is parallel to the Coriolis term when 
\u)\ increases (i.e. moving from the apogalacticon to p 6 rig cil etc t icon) and antiparallel when 
|a?| decreases (moving from the perigalacticon to the apogalacticon). Consequently, these 
latter two terms ere responsible of the initial de yiation of the tails from the r adial direction 
and of the formation of the S-shape profile (see lCapuzzo Dolcetta et al.l 120051 ) ). 



3.3. The tail direction along the orbit 

As a reliable quantitative definition of the tail orientation, we used the direction of the 
eigenvector corresponding to the greatest eigenvalue^ of the inertia tensor of a suitable por- 
tion of the cluster. More precisely, we considered t he stars contained i n an a nnulus centered 
at the cluster center-of-density (CD), as defined by lCasertano fc Hutl (119851 ) . and extending 



1 One of the three eigenvectors of the inertia tensor is parallel to the z-axis, while the other two lye in the 
(x, y) plane. Among these latter two, the eigenvector associated to the greatest eigenvalue gives the direction 
of maximum elongation of the stellar stream. 
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from « 30 pc up to a distance of 7-8 tidal radii (~ 170 pc). This region includes, indeed, 
the tail portion we call "inner" part of t he tails; thi s part is that usually detected in the 
observations. Note that in the study of lYim &: Led (120021 ) the direction of the tails was 
determined by eye, without any mathematically defined procedure. 



For each configuration of the system: 



• we evaluated the GC center-of-density; 

• we selected all the stars distributed between 30 pc and 170 pc from the cluster CD 
(excluding stars lying inside r < 30 pc, because this region corresponds to the spherical 
part of the system); 

• we evaluated the inertia tensor and the corresponding eigenvalues and eigenvectors of 
the selected stars; 

• we determined the eigenvector u corresponding to the maximum eigenvalue. 



This eigenvector is parallel to the direction of the two tails having opposite orientation. In 
the following, we always refer to the tail internal to the cluster orbit. 

The evolution of the tidal tail orientation is shown in Figs. El El 
These figures show the GC orbit and some kinematic quantities plotted as functions of time. 
From top to bottom we plot: the GC orbit (first panel), the distance r of the cluster from the 
Galaxy center (second panel), the magnitude of the angular velocity vector u> (third panel), 
the time derivative of |u?| (fourth panel). Finally, the fifth panel shows the evolution with 
time of two quantities: the angle 



_ -i / u ■ r GC . 
a = cos I - — - (7) 



l u ll r GC| 

between the tail direction and that of the galactic center and the angle 

p = cos 



between the tail direction and the GC velocity vector (the angles assume values in the 
[0°, 180° [ interval). 
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3.3.1. Tail-Galaxy center alignement 

In Figure [6], which refers to the orbit I, a is represented by two, barely distinguishable, 
solid curves, corresponding respectively to the GC with c = 1.02 and c = 0.8. The first result 
is that the time evolution of a is independent of the GC concentration (in the interval of c 
studied); this means that the morphology and orientation of the tails depend on GC orbit, 
rather than on GC internal parameters. 

Moreover, the tail is elongated toward the galactic center for most of the time: more in 
particular, if in Figj6]we exclude the regions around the pericenter, the average value of a 
is ~ 10°, i.e. the inner tail has only a slight deviation respect to the galactic center. 

When approaching the pericenter, instead, the Coriolis and the u> terms (in Eqj6]) are 
parallel and grow up. This induces a rapid alignement of the tails with the GC orbit. Beyond 
the pericenter, the Coriolis acceleration decreases and uj becomes antiparallel to the former. 
The net effect yields a realignment of the tails with the direction of the galactic center. 
Figure [3, which refers to orbit II, confirms the results found for the previous orbit. 

Finally, Figure [8] shows the case of the quasi-circular orbit (orbit III). Along this orbit, 
both the distance of the GC from the galactic center and the orbital angular velocity are 
nearly constant. This leads to a roughly constant a ~ 40°. This suggests that the change in 
time of the orbital angular velocity is what determine the tails orientation. 

3.3.2. Tail-orbit alignement 

The dashed curve shown in the bottom panel of Figures [6|71 and [HI represents the angle 
/3, previously defined. In the case of eccentric orbits, the first striking feature is that when 
approaching the pericenter the angle (3 decreases, reaches a minimum at the pericenter 
and then increases moving away from the pericenter. The maximum value of f3 is reached 
between the pericenter and the apocenter, depending on the shape of the GC orbit. However, 
in correspondence of the apocenter, f3 is ~ 90° in agreement with the value a ~ 0°. At the 
apocenter the tail is aligned along the radial direction and roughly perpendicular to the GC 
velocity. 

It is worth noting that, in all the cases considered here, the tidal tails deviate consider- 
ably by the GC velocity direction; only for very eccentric orbits and close to the pericenter, 
the angle f3 reaches a minimum value of ~ 17°. This indicates that the extrapolation of the 
cluster orbit from the elongation of this part of the tidal streams leads, in general, to predict 
GC paths with large errors. 



4. Conclusions 



This work is devoted to the study of the morphology and orientation of tidal tails 
surrounding globular clusters, in order to understand to what extent they trace the GC 
orbit and if some correlations exist between their orientation and the GC orbital phase. 
The main findings of our work can be resumed as follows: 



1. Tidal tails are good tracers of GC orbit only on the large scale, being the inner part of 
the streams never aligned with the GC path, unless the system moves on very eccentric 
orbit and close to the orbital pericenter. 

2. A strong correlation exists between the orientation of the inner tails and GC orbital 
position: they tend to be more elongated along the GC orbit when the cluster ap- 
proaches the pericenter, while, in turn, they tend to align toward the galactic center 
when the GC approaches the apocenter. 

3. We have shown that a key role in determining the tails morphology is played by the 
orbital angular velocity of the system and its variation with time. In the case of the 
least eccentric orbit, a nearly constant angular velocity determines an almost constant 
orientation of tails with respect to the cluster orbit and the galactic center direction. 
The amplitude of u> contributes, too, to the alignement of tails with the GC path, in 
particular just beyond the pericenter. 

4. The existing correlation between tails elongation and GC position along its orbit can 
be easily understood when referring to a non-inertial frame centered on the GC center, 
due to the role played by non inertial accelerations, as the Coriolis acceleration and 
that related to u) x r';. 



We want to outline that our findings are in good agreement with some o bservational 



result s. In particular we refer to the galactic globular clusters sample studied by iLeon et al. 



( 120001 ) , considering the 7 clusters for which the re liability level of the observed tidal streams 
is highest and orbital paramete rs are available (IDinescu et al.l Il999l ) . Six clusters of this 
subsample (see ILeon et al.ll2000l . Table 5) have tidal tails clearly elongated along the galactic 
center direction, and eccentricities greater than 0.62. This in agreement with our results 
that GCs, on eccentric orbits, have tidal tails elongated towards the galactic center for most 
of their path. The only exception is represented by NGC5139 (u Cen), which has tidal tails 
deviating from the galactocentric direction and extending towards the galactic plane; this 
situation does not contradict our results considering that NGC5139 is presently only 1.3 kpc 
above the galactic disk and is most probably undergoing a strong disk-shocking. 
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The simulations also repro duce the formati on of multiple tails around globular clus- 
ters, as observed for NGC288 (ILeon et al.l |2000| ) and for the galactic satellite Willman 1 
(jWillman et al.ll2006l ). Our simulations show that these features are expected for GCs on 
eccentric orbits near their apogalacticon, while they are absent in the case of GCs moving 
on less eccentric orbits. Actually, their formation can be explained by the behaviour of the 
Coriolis acceleration along the orbit: when its value is large, stars escaping from the cluster 
are pushed toward the orbital path. This is particularly evident around the pericenter; at 
the apocenter, where the Coriolis acceleration is weaker, stars tend to escape from the clus- 
ter more radially. At every pericenter passage, this effect produce a new tidal tail roughly 
aligned to the cluster path. 
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A. Equations of motion in the non-inertial reference frame 

Using a no-inertial reference frame (x', y', z') with the origin in the cluster center-of- 
mass, with the (x',y') plane coinciding with the cluster orbital plane and x' always pointing 
to the galactic center (so that this reference frame rotates with the GC angular velocity 
<jj respect to the inertial reference system (x,y,z)), the position of the i — th star in the 
galactocentric (x, y, z) reference frame can be expressed simply as: 

r» = r GC + r'i (Al) 

being r GG the position of the cluster center-of-mass in the inertial system. Rewriting the 
previous equation in terms of the (x', y') components (we omit the z' one for simplicity) we 
have: 

XiX + yiy = x GC x + y GC y + x\x' + y[y' (A2) 

being x,y and x',y' the unity vectors in the two frames, respectively. Derivating Eq JA2l with 
respect to time t, one obtains: 



XiX + yty = x GC x + y GC y + x\x' + y^y' + x^x' + y[y' = (A3) 
x GC x + y GC y + x'iX' + y\y' + u x (x-cc' + y\y') (A4) 
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being 

x' = u> x x' (A5) 
y' = uj x y' (A6) 

Derivating Eq JA3l with respect to time t and using Eqs. IA5HA6t one obtains: 

XiX + y\y = x GC x + y GC y + x'ix' + y' { y' + 2cj x (x'ix' + y^y') + (A7) 
cvx (u:x (x[x' + y'y)) + ux (x^x' + y[y') (A8) 

and so, finally, EqEJ 
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Fig. 1. — Orbits of a point mass with the initial position and velocity given in Table [3] in 
the MW potential (see sect. 12. ip . Top panel: orbit I; central panel: orbit II; bottom panel: 
orbit III. 



-17- 




Fig. 2. — Some snapshots of the GC moving along the orbit I (see Table [3] for the orbital 
parameters). One time unit corresponds to about 1 Myr. 



-18- 




Fig. 3. — Some snapshots of the GC moving along the orbit II (see Table |3] for the orbital 
parameters). One time unit corresponds to about 1 Myr. 
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Fig. 4. — Some snapshots of the GC moving along the orbit III (see Table |3] for the orbital 
parameters). One time unit corresponds to about 1 Myr. 



-20- 




Fig. 5. — Interpretation of the S-shape of the inner tidal tail around a globular cluster. 
The different terms in the rhs of the displayed equation (Eq. 7) are represented in the 
plot as arrows of different line styles. Note that the last term in the equation is here plotted 
antiparallel to the Coriolis term as it occurs when the GC moves from pericenter to apocenter 
(see text). The Galactic potential (included in the first and second term) is assumed, for 
simplicity, spherical. The cross in the lower part of the figure represents the Galaxy center. 
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Fig. 6. — Tails direction for the GC moving on orbit I. From top to the bottom: Plot of the 
GC orbit (some points along the orbit are marked with different symbols) ; the cross indicates 
the Galaxy center. 2: distance in kpc of the GC from the Galaxy center, as a function of 
time. The different symbols correspond to that in the previous panel. 3: GC orbital angular 
velocity in rad/Myr, as a function of time. 4: angular acceleration in rad/Myr 2 .5: solid 
curve: angle a formed by the inner part of the tails and the galactic center direction, vs 
time; dashed curve: angle /3 formed by the inner part of the tails and the cluster velocity, vs 
time. Both angles are in degrees. 



- 22 - 





8.- 



Same as FigEJ but for the GC moving on orbit III. 



Table 1. Parameters for the Galactic model (lAllen fc Santillan 



1991) 



Bulge 


M B 


1.41 x 1O 1O M 




b B 


387.3 pc 


Disk 


M D 


8.56 x 10 10 M Q 




a D 


5317.8 pc 




b D 


250.00 pc 


Halo 


M H 


10.7 x 1O 1O M 




a H 


12000 pc 



Table 2. GC structural parameters 

M GC [M Q ] r c [pc] r t [pc] c=log(r t /r c ) cr[km/s] 

4.7 xlO 5 2.2 23.5 1.03 10.5 

4.7 xlO 5 3.5 22.1 0.81 10.5 



Note. — Col. (1): cluster total mass. Col. (2-4): core 
radius, cutoff radius and concentration parameter of the 
King's model. Col. (5): central velocity dispersion. 



Table 3. GC initial orbital parameters . 
Orbit ID x\pc] y\pc] z[pc] ujkm/s] u y [km/s] ^[km/s] e 



I 


0. 


0. 


4000. 0. 


53.3 


0. 0.8 


II 


0. 


0. 


2500. 0. 


116.7 


0. 0.7 


III 


0. 


0. 


1500. 0. 


200.0 


0. 0.1 


Note. 


- Col. 


(1): 


orbit identification. 


Col. (2-7): 


cluster initial 



center-of-mass position and velocity components. Col (8): orbit eccen- 
tricity [Eq.(|SD] 



